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ABSTRACT 

Optically thin cooling gas at most temperatures above 30K will make conden- 
sations by pressure pushing material into cool dense regions. This works without 
gravity. Cooling condensations will flatten and become planar /similarity solu- 
tions. Most star formation may start from cooling condensations - with gravity 
only important in the later stages. The idea that some of the dark matter could 
be pristine white dwarfs that condensed slowly on to planetary sized seeds with- 
out firing nuclear reactions is found lacking. However, recent observations indi- 
cate fifty times more halo white dwarfs than have been previously acknowledged; 
enough to make the halo fraction observed as MACHOS. 

A cosmological census shows that only 1% of the mass of the Universe is of 
known constitution. 

Subject headings: White-dwarfs, Dark Matter, Cooling instabilities. 

1. Introduction - Russell in England 1902—1905 

After some time spent in attending lectures, Russell joined A. R. Hinks, a skilled 
observer, with whom he learned the techniques of accurate parallax determination. They 
used the Sheepshanks telescope, which cost 1 £3,183 in 1898. It was still at the Observatories 
for my early years in Cambridge, but was demolished in 1959. W. M. Smart, who used the 
telescope extensively, wryly remarked 

"It was a telescope of unusual design combining in a unique way the chief 
disadvantages of both the refractor and the reflector!" 

x To turn these prices into current dollars, multiply by about 100 
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In 1903 Hinks was promoted to Chief Assistant at a salary of £220 1 per annum and a house 
provided. Russell's Carnegie (post doctoral) grant was £200 1 per annum. 

It is a tribute to the persistence of Hinks and Russell that they eventually determined 
50 parallaxes, many of which stand up well to comparison with modern determinations. A 
not atypical example is the parallax of the well known high velocity subdwarf Groombridge 
1830 for which Russell gives 0.100 ± 0.029" as compared to a modern value of 0.107". 
Thus, many of the stars in Russell's first HR diagram were put there on the basis of their 
Cambridge parallaxes. 

2. Census of the Universe 

The prediction that the angular wavelength of the first peak in the angular power 
spectrum of the cosmic microwave background would decide whether the Universe was 
closed or open was made only 14 years ago (Efstathiou & Bond 1986; Bond & Efstathiou 
1987) [see also Doroshkevich et al. (1978)]. The first results of TOCO, Boomerang and 
Maxima are already with us. By contrast the prediction that many large galaxies including 
our own, M31, M32, M81, M82, M87, etc, contain giant-black-hole remnants of dead 
quasars in their nuclei (Lynden-Bell 1969) took over 26 years before the first definitive case 
was found (Miyoshi et al. 1995). Although there were indicative results earlier (Young et 
al. 1978; Sargent et al. 1978; Penston et al. 1981; Pounds et al. 1986; Dressier & Richstone 
1988; Kormendy 1988) and many convincing ones since (Eckart & Genzel (1996); Genzel et 
al. (1997); Gebhart et al. (2000)). An interesting account of twentieth century developments 
is given in Shields (1999) review, but it can be argued that the first prediction was 211 
years earlier (Michell 1784). 

Astronomy is now running so fast that almost everything is as ephemeral as a news 
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report, nevertheless each entry in Table 1 is attested by at least two independent lines of 
evidence. For example, the cold dark matter fraction comes both from estimates of the 
baryon fraction in large X-ray emitting clusters of galaxies (White et al. 1993) and from 
estimates of the relative height of the peaks in the Maxima results. Likewise the unknown 
expansion energy estimate comes both from the difference between the closure density 
(TOCO, Boomerang and Maxima) and the total matter density of Maxima, and from the 
rate of acceleration of the Universe measured using distant supernovae. 

However, if in a thought experiment we dream of reversing the Hubble expansion - 
but nevertheless leaving time running forward - then the black body radiation from the 
sky would get hotter and hotter. The white dwarfs, neutron stars and black holes would 
all survive to temperatures of more than 10 9 K and there is then no time to destroy them 
before the big crunch. This thought experiment shows us that any compact bodies which 
might already be present at the time of Helium creation 10 seconds after the Big Bang 
would survive to the present day. If so, their baryons would not be counted in the first two 
entries but would contribute to the collisionless matter of the third entry. The two entries 
for baryons disagree but the second estimate is based on the results of Boomerang and 
Maxima which don't agree well either. Few believe the number of baryons goes up between 
10 sec and 1/3 Myr and the first entry is considerably more secure. However, even that 
gives us a missing baryon problem because the numbers of baryons in stars, galaxies 
and intra-cluster gas is much less than that given by the first entry. Much may lie as yet 
undetected as 5.10 5 K gas in small groups of galaxies like ours. Detection of the dispersion 
measure to pulsars in the globular clusters of the Fornax dwarf would test this. 

The wonder of astronomy lies as much in what is unknown as in the beauty of the 1% 
that is known. However, when discussing these unknowns, one should heed the wisdom of 
the past, so when lecturing about cold dark matter, I always quote 
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TIME FRACTION = p/p c Q 

H' 1 = lO^/i-Vrs p c = 3H 2 /(8nG) h = 0.6 

lOsec Baryons (Atoms) 1.9h' 2 % 5% 

±0.1 

1/3 Myr Baryons (Atoms) 3.2h~ 2 % 9% 

±0.6 

Dark Matter ± Baryons 15/i~ 2 % 40% 

±5 

10 10 yr Radiation 0.007% 

Neutrinos 1% ? 
Unknown Expansion 

Energy f2 A (Supernovae) ~ 60% 

1/3 Myr Closure 1.03 ±0.1 100% 



Table 1: Universal constituents and effective times of measurements. 
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"If a thousand men believe a foolish thing it is still a foolish thing." 

A new way of looking for invisible heavy halo objects via gravitational lensing events towards 
the Magellanic Clouds was first proposed by Petrou (1981). Independent calculations were 
made by Paczyhski (1986), and it was he who stimulated Alcock et al. (1996) who first 
convinced the world that it was a practical possibility and led the Macho team to find such 
events (Alcock et al. 2000). 

The Macho team found over twenty events toward the LMC which could account for 
between 10% and 20% of the Galaxy's dark halo, although almost half of the events might 
be self-lensing within the LMC. Whereas they found too few events to account for their 
halo model, both they and the Eros and Ogle teams found more than the expected number 
of lensing events towards the Galaxy's bulge. It would be not unnatural to imagine that 
the models of the Galaxy's mass distribution were seriously wrong with too little in the disc 
and too much in the halo, but that is not how the experts put it. 

While pondering the nature of dark matter and the Macho results on Feb 15th 2000, I 
realised that there was a possibility that I had never heard discussed, 

"Could white dwarfs be made directly from primordial material without ever 
forming stars or burning hydrogen? Could such pristine white dwarfs be the 
halo objects found by Macho?" 

As we shall see, such questions had occurred earlier to Salpeter (1992) and Hansen (1999). 
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3. Pristine and Halo White Dwarfs 

Discussions of star formation often start with a large body contracting from a diffuse 
cloud. The equilibrium radius of such a body is determined by the Virial theorem. 

r=-iv = -*, 

where T the kinetic energy is ^NkT and V the potential energy is —GM 2 /(2R). The 
radiative loss from the surface causes the energy to get more negative so the radius R 
decreases and the internal temperature T increases until the nuclear reactions fire and the 
body becomes a star. 

The above conclusion can be circumvented if the kinetic energy is not due to the star's 
temperature but is determined by the zero point energy of confined electrons required by 
the uncertainty principle. This occurs in solids and liquids. 

If I were to freeze the A.A.S. President to absolute zero, she might become glassy eyed 
and she would shrink by an inch or so, but not substantially. This is because the electrons 
in her atoms owe their kinetic energy neither to their temperature nor to their angular 
momentum (since there is none in e.g., the ground state of hydrogen) but to their zero 
point degeneracy energy. They are of course held together by electricity. 

Asteroids, planets and white dwarfs all owe their sizes to this zero point energy of 
electrons. In fact due to the work by Fowler (1926), Stoner (1929), Anderson (1929) and 
Chandrasekhar (1931) and by Salpeter (1967) in the 1960's, the mass radius relationship 
for such 'cold' bodies is one of the best understood parts of astrophysics. 

In outline there are three terms: 

• The zero point energy of the electrons which may or may not be relativistic. 

• Their electrical binding to their associated atomic nuclei. 
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• The gravitational potential energy. 

In small bodies such as asteroids, the main balance is between the zero point energy and 
the electrical attraction of the nuclei. Gravity merely serves to keep the parts together 
so as the bodies are piled together the mean density remains much the same. In spite of 
the increased pressure at the centre this equal density law still holds up to masses close 
to Saturn's. About there the gravitational binding starts to influence the total binding 
energy which is still mainly electrical. As a result the mean density starts to increase. At a 
mass a little over the mass of Jupiter gravitational and electrical binding become equal and 
the weight of each extra mass added causes a contraction as great as its volume. Beyond 
that point the bodies get smaller as the mass increases. When the zero point energy of 
each electron is substantially less than its rest mass, we get the brown dwarf-white dwarf 
sequence with R oc M^ 1 / 3 . Once the electrons become relativistic the pressure generated 
becomes as soft as gravity and the radius decreases very rapidly as the Chandrasekhar 
mass is approached (Anderson 1929; Stoner 1930; Chandrasekhar 1931). It is possible to 
get a single mathematical formula for the whole mass-radius relationship from asteroids to 
Chandrasekhar's limit (Lynden-Bell & O'Dwyer 2001). 



where M p is the mass of the planet of maximum radius and p is the density of the 
solid at zero pressure. I is unity except close to the Chandrasekhar limit M c h where 



The mass of the planet of maximum radius is set by the equality of gravitational and 
electrical attractions. 



The Chandrasekhar mass is approximately (hc/e 2 ) 3 ^ 2 times greater than this, showing that 
the fine structure constant plays a vital role. 




/ = J/[l + (3/4)(l - J 1 ' 2 )} with J = y/l - {M/M ch y/\ 
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In the theory of gravitational fragmentation the mass of the non-fragmenting fragment 
is given by (Low & Lynden-Bell 1976) 



M 



3 5 • 5 4 / hc\ ( he \ 10 / m e 



1/7 

m H ~ M /6O , 



2 8 7r 5 \e 2 J \Gm 2 H ) \ m -H\ 

which may determine the minimum mass of bodies so formed, but as yet there is little 
evidence that nature pays too much attention! See in particular the fascinating a Orionis 
Cluster (Zapatero-Osorio et al. 2000). 

Figure 1 shows two mass-radius relationships; that on the left is appropriate for rock or 
for white-dwarfs with (// = 2) two baryon masses per electron. That on the right (/i = 8/7) 
is appropriate for a pristine mixture of 75% Hydrogen and 25% Helium by mass, i.e., for 
Hydrogen planets, brown dwarfs and - if such exist - pristine white dwarfs. The latter 
have radii R oc /i~ 5//3 M _1 that is 2.5 times larger than equal mass normal Helium, Carbon 
or Oxygen white dwarfs. At the same surface temperature we would expect them to be 
two magnitudes brighter. Since they have twelve times as many nuclei as Carbon white 
dwarfs their specific heats will be greater and they may not cool faster in spite of this extra 
emission. 

When first interested in pristine white dwarfs I consulted early works by Parenago 
(1946) and by Eggen & Greenstein (1965) which showed two sequences separated by nearly 
two magnitudes. It was rather ironical to find that after Eggen consulted a theorist his 
work on the Hyades white dwarfs restricted itself to candidates within one magnitude of the 
expected sequence and would therefore have missed a second sequence. The idea was that 
in low-density regions such as cooling flows, galaxy haloes and dwarf spheroidals, stars form 
very gradually by accretion on to planetary sized seeds. If the growth was so slow that the 
accretion energy was radiated before the material was buried then the body would grow in 
mass resting on its zero point energy and never getting hot enough to fire its nuclear fuel. 
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When L-B raised the question of making pristine white dwarfs with C. A. Tout, his 
immediate reaction was that they would blow up like novae, but a day later he had found 
that novae only explode when the temperature of the degenerate hydrogen layer exceeds 
10 6 K. If the accretion rate were so small that the temperature never reached a million then 
pristine white dwarfs might be made. 

To see if any had been observed we next consulted Mike Irwin. He said "it is possible 
that we have just what you are looking for" and sent us to read the growing literature on 
halo white dwarfs. 

When Ibata et al. (1999) repeated the Hubble deep field three years later they found 
five objects at 29 th magnitude which appeared to move. One may be a distant supernova in 
the side of a pre-existing image, two of the other images have unexpected asymmetries, but 
two looked like single stars 3 that move 20 milliarc sec/yr. 

If such a star were brought 100 times closer it would be at 19 th magnitude and move 
2 arc sec/yr. Thus Irwin planned a wide-field Schmidt Survey to find such rapidly moving 
objects from existing plate material with plates taken less than ten years apart so that the 
faint stars would not move so far that they were 'lost' among others. 

After searching over 700 square degrees Ibata et al. (2000) found 18 possible objects 
of which the brightest two were investigated spectroscopically. Both showed the strong 
infra-red blanketing predicted as the hallmark of molecular hydrogen formation in cool 
white dwarf atmospheres (Mould & Liebert (1978); Hansen (1998)). Both had motions 

3 Since this lecture was given, a third epoch HST frame has shown that neither of these 
remaining objects has a uniform motion, but the objects found by Schmidt surveys are real 
and more candidates have since been found (see Oppenheimer et al. (2001) and the criticisms 
of Reid et al. (2001)). 
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appropriate for halo stars 'left behind' by the rotation of the Galaxy. Earlier Hodgkin et al. 
(2000) had found another white dwarf showing the strong predicted infra-red blanketing. 
To be so fast moving and so dim the objects cannot be main-sequence stars without leaving 
the Galaxy. Making them old white dwarfs means they are quite close, about 30 parsecs. 
This means such objects are common; 50 times more common than earlier estimates of halo 
white dwarfs, but not all estimates (Chabrier et al. (1996); Chabrier (1999)) showed that if 
the halo mass function peaked a little above 1M then there could be many cool old white 
dwarfs. At such a density they would account for up to 10% of the dark halo. Others have 
objected that the population of heavier stars born with them would have produced more 
metals than are observed, but that would be avoided if these were pristine white dwarfs. 
At this juncture it seemed that the dark curtain over the nature of dark matter had let 
through a chink of light. L-B started to study how stars might form at low densities in 
cooling flows but, in the meantime, Tout was in real difficulties. Pristine white dwarfs over 
about 0.2M Q always fired their hydrogen unless their accretion rates were only 1 M Q /10 13 
years - too slow to be made in this Universe. Furthermore this agreed with the estimates 
made by Salpeter (1992) on the basis of the accretion models of Lenzuni et al. (1992). Just 
as nuclear physicists predict an island of almost stable ultra-heavy nuclei that are very hard 
to make, so pristine white dwarfs between O.2M and the cold pycnonuclear reaction limit 
of 1.1 M Q are a theoretical possibility, but no-one knows how to make them! 

This makes yet more interesting the origin of the large number of halo white dwarfs 
now being found. (See also Harris et al. (2001)). 

4. Cooling Instability 

Although the investigations of sections 4-7 were stimulated by the consideration of 
pristine white dwarfs, they are independent and develop the non-linear theory of cooling 
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instabilities which are probably important precursors of galaxy formation, star-cluster 
formation and normal star formation. 

As stars are held together by gravity it is natural to believe that gravity is the primary 
driver in their condensation, but this may be wrong. At a density of 10~ 22 g cm -3 a solar 
mass of gas has an escape velocity of only 0.1 km/s whereas an ionised medium of that 
density will have a sound speed of about 10 km/s. At low densities pressure can be much 
more important than gravity and the pressure is often dictated by the cooling as the gas 
radiates to infinity (Fabian 1994). 

Anyone who looks at photographs knows that the idea of a uniform interstellar medium 
is a theoretical abstraction far from truth. With such variations from place to place the 
cooling rate per gram will vary. This will induce pressure differences that will push the 
interstellar gas to the low pressure regions. The cooling is commonly such that regions of 
lower pressure have higher density. Thus condensations can form due to pressure alone, 
quite independently of any gravity. This is how condensation begins, the gravity comes in 
later. 

We shall consider the simplest possible case of an optically thin gas at pressure p 
cooling by radiation to infinity. Suppose there are regions of different densities with 
anti-correlated temperatures so that dynamical equilibrium is maintained with the pressure 
spatially uniform. The second law of thermodynamics gives the entropy decrease consequent 
on the cooling p 2 A(T) per unit volume; so for unit mass 



T Ds/Dt 



pA(T) . 



(1) 



For a perfect gas, s = km 1 ln(T 3 / 2 /p) + constant = km 1 ln(T 5//2 /p) + constant where m is 
the molecular mass so for two regions of gas at a common pressure p = m~ 1 kTp, 



(5/2)fc 2 m~V 1 J D/ J Dtln(T 1 /T 2 ) = T~ 2 A(T 2 ) - Tf 2 A(T X ) . 



(2) 
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Hence as the gas cools the temperature ratio will increase whenever T~ 2 A(T) is a decreasing 
function of temperature. This is the cooling instability criterion discovered by Fall & Rees 
(1985). As they pointed out it differs from the T _1 A(T) decreasing criterion for thermal 
instability found in Field's masterly paper (Field 1965) and originally due to Weymann 
(1960). The reason for the difference is that Weyman and Field studied small departures 
from equilibrium conditions in which heating and cooling balanced in the mean state while 
the cooling instability criterion holds when everything is cooling on a timescale that is slow 
enough that the two regions share a common pressure (which may vary with time). 

Figure 2, adapted from Dalgarno & McCray (1972), shows that for 30K < T < 9000K 
and T > 10 4 K the whole range in which gas is ionised, T~ 2 A(T) decreases. We deduce 
that in the absence of heating, cooling instabilities will occur and cooling condensations 
will develop. However, we have neglected two effects that work to suppress the instability. 
At small scales thermal conductivity will tend to suppress temperature differences and 
thus depress the instability by heating the denser regions. But, once the instability is kick 
started the thermal conductivity falls as the temperature falls so limiting the effect. More 
importantly Balbus & Soker (1989) have shown that, for infinitesimal perturbations, there 
are no thermal instabilities when the gas is in stratified equilibrium in a gravitational field. 
How can this be? Surely if we make the gravity field weaker and weaker we must arrive 
at the uniform medium instability criterion once the gravitationally induced gradients 
are weak enough! Here it becomes necessary to contrast mathematically exact results for 
infinitesimal perturbations with physical results for finite perturbations (Lufkin et al. 1995). 
An Astrophysicist might well consider a perturbation with Sp/p of 1/10 or even 1/3 to 
be quite small, but a Mathematician will consider only the limit as Sp/p tends to zero so 
that all its gradients are small compared with those that he has in his unperturbed state. 
For weak gravity fields this is a severe restriction on perturbation amplitudes especially 
at short wavelengths. We find in the next section that new density maxima caused by 
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the perturbation are all important. In an infinitesimal perturbation new density maxima 
only occur in a uniform medium. However small the gradient in the unperturbed state, 
infinitesimal perturbations make no new density maxima. For systems with a gravity field 
our results will only apply to finite amplitude perturbations which make new maxima. 



5. Non-Linear Cooling Condensations 

We take the simplest possible case of no external gravity, no self-gravity, no thermal 
conductivity and slow cooling so that the pressure is spatially constant. We consider this 
not because it is the most realistic case but because it gives insight and insight rather 
than complicated exactness is what science is about (Eddington 1926). Once the insight 
is obtained more complications are more readily understood. For a plasma cooling by 
free-free emission the cooling rate per gram can be written Kpc s where c s the sound speed 
is proportional to T 1 / 2 . The coefficient K has the dimensions [M _1 L 4 T~ 2 ] which is one 
power of L more than the dimensions of G. As an aside it is interesting to ask how long this 
length K/G is, since both free-free emission and gravity are inevitably involved in galaxy 
formation. In fact in fundamental dimensionless constants and the classical electron radius, 

3 \ 3 m p J \nc j \Gm e m p J m e c z 

It is this large distance or rather 1/3.8 times it, that gives the 73 kpc radius from which 
galaxies fall together because their cooling rate is faster than their collapse time. We wish to 
discuss more general cooling laws but, although it is not hard to be more general, simplicity 
suggests a general power law oc pT 2 ~ a so the above case corresponds to a = 3/2. With our 
general power law equations (1) and (2) together with the perfect gas law take the form 

5/2 D/Dtln (p 3/5 /p) = -K (p/p 3/5 ) a p 1 ' 2 ^ 5 , (3) 
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so 

D/Dt (p 3/5 /p) a = -(2/5)aKp 1 ' 2a / 5 . (4) 

Writing p for the initial value of p on the fluid element considered and momentarily 
considering the special case with p independent of time we find 

(pS/B/ p )« = ( p 3/5 /po) * _ (2/5)aKp 1 - 2a /H 

and thence 

p/ Po = [l-(2/5)aKp 1 - a p^t]- 1/a (5) 

However (4) may also be solved when p depends on time. We introduce a new weighted 
time r defined by 

t= J(p/Poy- 2a/5 dt, 

where p = p(t) and p is its initial value. 

Then in place of (5) we have the general solution 

p/po = ( P /Po) 3/5 [1 - {r/r c )]- 1/a , (6) 

where r c = 5/(2o;_ft'po~ a Po) i s the 'time' at which the density of our fluid element become 
formally infinite. With the Fall-Rees criterion satisfied a > so this time is shorter when 
the initial density po is greater. Thus nucleation occurs around maxima in the initial 
density distribution. We rewrite (6) in terms of r m the collapse time for the initial density 
maximum p m . Multiplying (6) by po/p m we have 

P/Pm = (p/Pof /5 [(po/ Pm )- a ~ r/r m ] ~ 1/a . (7) 

At time r the current maximum density of the material whose initial density was p m we 
call Pcij) or p c for short. For the maximum density we put p = p m in (7) and have 

Pc/p m = (p/po) V5 \l-r/r m r 1/a (8) 
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for t < r m . Of course for r > r m the maximum density is infinite but we find it convenient 
to define p c as a characteristic density given at all times by (8) even when r > r m . This 
characteristic density is the maximum density at all times prior to r = r m and increases. 
After r = r m it becomes smaller again and represents the density at a self-similar point in 
the ensuing accretion flow. Dividing (7) by (8) we obtain 

-l/a 



P_ 

Pc 



(Po/Prn) °-\± 1 



(9) 



|1 -r/r, 

where the plus sign is to be taken for r < r m , before the centre collapses and the minus sign 
afterward. As r — > r m the denominator becomes small, so (9) is then very sensitive to small 
differences between p an d p m - This means that the density profile near the collapse time 
r m depends solely on the initial profile of p close to p m . We now show how this leads to 
similarity solutions. 

Let M(> p ) be initially the total mass, centred on one maximum of density p m , which 
has at each point a density greater than some chosen value po- For a spherical density 
distribution around the maximum we have 

dM = A7ir 2 p dr = (dM/dp )dp . 

For a cylindrical distribution we replace 4nr 2 p dr by 2nRp dR while for a planar 
distribution we would replace it by 2p dz and measure mass per unit area. Thus 

(4/3)vrr 3 
irR 2 



2UI 



/ {p' Q )-\dM/dp' )dp' Q = [po(M')]- 1 dM' , (10) 

Jpo Jo 



where at the last equality we have used not M(p ) but its inverse function po(M) which 
gives the density at the edge of the mass M of higher density. Now according to our cooling 
law higher densities get denser faster, so the density ordering of fluid elements remains the 
same. If a given initial density p evolves into a density p after 'time' r then the function 
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Po(M) will evolve into p(M) with the same M. Thus the radius of the mass M at time r 
will be given by analogy with (10) as 

(4/3)7rr 3 



nR 2 
2\z\ 



M 



l/p{M')dM' 



Now near any smooth quadratic maximum the initial density takes the form 

Po = Pm/(1 + r 2 /a 2 ) = p m (l + ar 2 /a 2 )~ 1/a , 
[where r 2 should be replaced by R 2 or z 2 for the 2D or ID cases]. 
Hence 

M = (4/3)7rr 3 p m [l + 0(r 2 /a 2 )} , 
so eliminating r in favour of M 

' p m {1 + a[3M/(4ira 3 p m )] 2 / ? >y 1/a spherical 
Po(M) = 1 p m {l + aM/(na 2 p m )y 1/a cylindrical 
p m {l + a[M/(2ap m )] 2 y 1/a planar 
Inserting this into (9) we find 

f p c (m 2 / 3 ± l) 1 /" 
p (M) = I Pc (m ± i)V« 
k p c (m 2 ± lf/ a 



> = Pc{r)p*{m) 



where 



3a 3/2 M/(47ra 3 p m |l - r/r m | 3 / 2 ) spherical 
in = < aM/(7ia 2 p m \l — r/r m |) cylindrical 
w a 1/2 M/(2ap m |l -r/r™] 1 / 2 ) planar 



and 



p* 



(m 2//3 ± l) -1 /" spherical 
(m ± l) -1 /" cylindrical 
(m 2 ± l)~ 1/a planar 



(11) 



(12) 



(13) 



(14) 



(15) 



(16) 
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To find r we insert (14) into (11) and deduce 

(4/3)7rr 3 p c = / (m 2 / 3 ± l)^ a dM (17) 

Jo,Mc 

where the lower limit is zero before collapse and the collapsed mass M c afterward, as that 
collapsed mass does not contribute to the radius. Using the dimensional pieces of (15) and 
(8) to define a characteristic radius r c (t) we find r = r c {t)r\{m) where 



r 3 



/ (m 2/3 ± l) 1/a dm , 
Jo,i 



similarly 



Rl = 


J a+l 1 

I aTi( 

— 


■ = N/ 


' z c and 




= a\l 


Rc{r) 


= a\l 




= a 1 



l~l) l+1/a T<T m 

f™(m 2 ± l)V«d m 

-T/T^V2+l/3a a -l/2 (po/p) l/5 

-r/r m |V2 + V2a a -i/2 (po/p) 3/io 
■r/r m | 1 / a + 1 /« a -i/2( po / p )S/5 



(19) 



at least when p is constant. 

z c (t) < R c {j) < r c (r) as r approaches r m , and (19) always holds near r = r m since 
p(r) is not sensitive to r m . Thus the planar collapse solution collapses fastest, just like the 
gravitational case. This strongly suggests that the spherical solutions are unstable to a 
flattening instability and we shall see that this is indeed true. 

Equations (16) and (18) together constitute parametric equations for the density 
profile p*(^*) in terms of the parameter m. From (14) it is already clear that the solution 
is self-similar both before and after the moment when the density first becomes singular. 
The spherical solutions with a = 1 are plotted in Figure 3. This self-similarity rests on the 
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approximate form of po close to its maximum p m , so it only holds for times close to that 
moment. One may find p*(r*), or rather its inverse function r^p*), explicitly for the special 
cases when a = 2/n, with n an integer, but we keep to the general case and look at the 
asymptotic forms for small and for large r*. 



Pre-collapse 



For r* small r 3 = m + (3/5)a _1 m 5 / 3 so 

m = — (3/5)o;~ 1 r^ and from (16) 

p* ~ (1 + r 2 )~ 1/o r* small 



(20) 



for r* > 1 



'* 3a+2'" 



3a 2/(3a) 



[m + 



2_t3a m l/3 

2+a 



...] 



m — T] ^rf 1 



1 " 3T^"V +2, ' /3r *~ 2, ' + °( r * 



(21) 



where rj = 3a/(3o; + 2) 

so for large r*, (16) gives p, = (1 + ^-Vs^-i/a K r -e/(3a+2) 

Thus at the moment when the central density becomes singular we find a power law density 
profile with index -6/(3a + 2). That is -12/13, -6/5 and -12/7 for a = 3/2, 1 and 1/2 
respectively. 

The r~ 12 / 7 profile is the same as that for Penston's cold free-fall solution for a self gravitating 
spherical cloud, (Penston 1969). 



Post collapse 



For times later than r m the central point mass M c grows 
M c = [(r/r m ) - if 2 (4/3)7ra 3 p m «- 3 / 2 



(22) 



which comes directly from expression (15) with m — 1. Of course (22) only holds when r 
is close to r m when approximation (12) is valid. However, it is easy enough to determine 
M c (r) directly from the initial conditions. Let M[p ] be the mass within the po density 
contour in the initial conditions and suppose it all has density > p for all p above some 
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minimum. Then from our expression for r c above (7) the mass that has collapsed by time r 
is the mass that had initial density greater than [(2/5)aKpl~ a r] 



-l/a 



i.e., 



M r = M 



{(2/5)aKpl- a r) 



-l/a 



(23) 



Setting m = 1 + rh in (18), the post collapse similarity solution for times a little greater 
than r m is (see Figure 3c) 

r 3 = (2/3) 1/Q [a/(a + l)]m 1+1/Q + ... 

p- a = (2/3)m + ... (24) 
hence p* = {(3/2) [a/(a + 1)] r^} 1/{a+1) + . . . small r* 
i.e., a power law of index —6/5, —3/2, or — 2 for a = 3/2, 1, 1/2. At r > r c (t) these 
power laws change over to the power — 6 /(3a + 2) derived earlier. 



6. Exactly Self-Similar Solutions 

We derived these solutions as approximations, only good for times close to r m for 
regions not too far from the maximum density. However if, instead of postulating a general 
initial profile and then approximating it near maximum by (12), we had taken an initial 
profile defined by (18) and (16) then our solution would be exactly self-similar at all times 
and all positions. Perhaps the simplest of these is that given by Newton's law of cooling 
proportional to temperature. Then (18) integrates and yields for a — 1 

rl = m + (3/5)m 5/3 r < r m 

(3/5) (m 5 / 3 - 1) -m + 1 r > r m 

while (16) yields = l/(m 2 / 3 ± 1) 

Asymptotically p* — > r* 6 ^ 5 for all large r* so to keep the pressure constant T — > r*^ 5 . 
Similarly the 'exact' solution for the cylindrical case gives p explicitly 

P = p c (t)/[U+(l + a- 1 )Rl] l/{a+1) 
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where U = 1 for r < r m and otherwise; p c (r) is still given by (8) in all three geometries. 
For « = lwe have for the planar case 
p, = l/(m 2 ± 1) 
and z* = \z\lz r = 



(l/3)m 3 + m r < r, 

(l/3)m 3 -m + 2/3 r > r, 



We have obtained all these solutions under the assumption that the cooling is so 
slow that we may take the pressure to be almost constant - i.e., the system evolves 
quasi-statically through a sequence of equilibria. Nevertheless if it evolves in finite time 
there must be small motions and very small accelerations to drive these small motions. To 
find these motions and the small pressure gradients that drive them we use the equations 
of fluid mechanics. By mass conservation 

D\np/Dt = — div u = —r~ 2 d/dr(r 2 u) 

where u is the radial component of u. We know D\ap/Dt from (3) so using (18) 



u = -(l/5)(p/p)r - {2/lh)Kp l - a p a c r c r- 2 



Jo,i 



When a = 1 the final integral simplifies to m for r < and m — 1 for r > 0. Such 
solutions are plotted for p = as Figure 4. 

In 2D 

(m+ l) 1/a - 1; r < 
( m _ 1)V« ; r > . 



u = -(3/10) (p/p)r- (l/5)Kp 1 - a p^R c r- 1 a 



In the planar case with a = 1 

u = -(3/5){p/p)z - (2/h)Kp l - a p a c r c m 
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Notice that similarity only extends to the velocity fields when either p is zero or 
p oc p 2 ~ a p". From the above Du/Dt, and thence the small pressure gradient can be 
calculated. Spherical self-similar solutions for the outer parts of cooling flows have been 
found numerically by Bertschinger (1989). His solutions are valid through the sonic region. 

7. Flattening into Sheets 

We saw in section 5 under (19) that the collapse time depended on the dimension and 
the planar collapse was fastest. We were thus led to believe that density maxima would 
flatten and the cooling instabilities would lead to curtain-like sheets of high density. To 
check that this is indeed the case, Marcus Briiggen, ran some 3D and 2D simulations. In 
these the full equations of fluid dynamics were integrated and the cooling was switched 
on over a few sound-crossing times to avoid transients. While the pressure was initially 
constant in space it was allowed to vary in space and time. While the cooling was slow it 
was not so slow that pressure remained constant and indeed it fell drastically in the cooling 
sheets that developed. Figures 5a and 5b which contour the density distribution initially, 
and after cooling, show clearly the development of flat sheets of high density just like those 
studied in planar geometry by Burkert & Lin (2000). The high density and low sound speed 
in these sheets make them ideal for gravitational instability and fragmentation into clusters 
of stars. 

8. Conclusions 

While some sort of census of the contents of the Universe is now possible, the true 
nature of the vast majority of it is unknown. Astronomy is still young. 

Many more halo white dwarfs are now being found and this is likely to change our 
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concept of what the galaxy is like. They probably make a significant contribution to 
the dark matter (5% - 10%), however there remain serious doubts on the reality of this 
interpretation. This is an exciting and controversial subject that deserves our attention. 

Pristine white dwarf configurations exist for objects that never burn hydrogen but 
no-one has devised a way of getting bodies to that state so the large numbers of halo white 
dwarfs now being found are probably normal burnt out stellar remnants. 

Cooling instabilities may be the method by which density enhancements start in many 
different areas of astronomy; gravity may be only important in the final stages of star 
formation. 
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Fig. 1. — Figure 1. The heavy line on the left gives the mass-radius relation for cold bodies 
with p = 3.65 g cm" 3 and /i — 2, suitable for rocky planets and normal white dwarfs. 
The heavy line above gives the black hole radius. The lighter line is for a mixture with 
75% Hydrogen and 25% Helium, fi = 8/7 suitable for gas-rich planets, brown dwarfs and 
pristine white dwarfs. This sequence ends at 1.1M due to pressure induced, (pycno-nuclear) 
reactions. Stable bodies could exist between 0.2 and 1.1M© but no plausible formation path 
has been found. The dashed lines extrapolate the non-relativistic formula for degeneracy 
into the region of its invalidity. The neutron star line is parallel to these but displaced from 
the pure hydrogen line by the factor (m e /m n ) in radius. 
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Fig. 2. — Figure 2. The cooling function A(T) of Dalgarno &McCray with superimposed 
(dashed) the graph of 10 8 A(T)/T 2 . There is cooling instability causing condensations in 
an otherwise homogeneous region wherever the latter graph is falling, i.e., everywhere with 
T > 30K except the region near 10 A K where hydrogen is partially ionised. Cooling instability 
should be prevalent throughout astronomy. 
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Fig. 3. — Figure 3a. The density profile for the spherical self-similar solution with a — 1. 
Figure 3b. The log density - log radius profile at times — = —1, —10, —100 showing how 
the scale of the profile shrinks as r = is approached. Figure 3c. The same as 3b but after 
the centre collapses. At r = 1 there is a growing central mass and the radius at which the 
power law changes from —3/2 to —6/5 increases like (r/r m — l) 5 / 6 again for a — 1. 
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Fig. 4. — Figure 4a. The inward radial velocity profile prior to collapse of the core. Figure 
4b. The same after the point mass develops u oc r~ 1//2 for small r and u oc r -1 / 5 for large r. 
The transition point moves outward. 
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